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CORRELATION OF FINITE ELEMENT FREE VIBRATION 
PREDICTIONS USING RANDOM VIBRATION 
TEST DATA 

JEFFREY A . CHAMBERS 
ABSTRACT 

Finite element analysis is regularly used during the 
engineering cycle of mechanical systems to predict the 
response to static, thermal, and dynamic loads. The finite 
element model (FEM) used to represent the system is often 
correlated with physical test results to determine the 
validity of analytical results provided. Results from dynamic 
testing provide one means for performing this correlation. 
One of the most common methods of measuring accuracy is by 
classical modal testing, whereby vibratory mode shapes are 
compared to mode shapes provided by finite element analysis. 
The degree of correlation between the test and analytical mode 
shapes can be shown mathematically using the cross 
orthogonality check. 

A great deal of time and effort can be exhausted in 
generating the set of test acguired mode shapes needed for the 
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cross orthogonality check. In roost situations response data 
from vibration tests are digitally processed to generate the 
mode shapes from a combination of modal parameters, forcing 
functions, and recorded response data. 

An alternate method is proposed in which the same 
correlation of analytical and test acquired mode shapes can be 
achieved without conducting the modal survey. Instead a 
procedure is detailed in which a minimum of test information, 
specifically the acceleration response data from a random 
vibration test, is used to generate a set of equivalent local 
accelerations to be applied to the reduced analytical model at 
discrete points corresponding to the test measurement 
locations. The static solution of the analytical model then 
produces a set of deformations that once normalized can be 
used to represent the test acquired mode shapes in the cross 
orthogonality relation. 

The method proposed has been shown to provide accurate 
results for both a simple analytical model as well as a 
complex space flight structure. 
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CHAPTER I 


INTRODUCTION 


Ltl Background 

During the engineering development cycle of 
structural/mechanical systems, finite element models are often 
used to provide insight into the static and dynamic response 
characteristics of the system. Given the broad assumptions 
made during the finite element modeling process, there is a 
need to verify the accuracy of the modeling technigues 
employed. The validity of these math models is usually 
substantiated through static or dynamic tests conducted on the 
engineered system. Results taken from these physical tests 
are compared with analytical results provided by the finite 
element model (FEM) under the same boundary and loading 
conditions. 

Probably one of the most common methods used to validate 
a finite element model is through the use of modal testing and 
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correlation [1], With this technique, the physical test 
specimen is fixtured to represent its working environment, or 
may be suspended in such a manner as to simulate free 
restraint condition, and is subjected to a known forced 
excitation (or simultaneous excitations). The response of the 
system to this excitation is measured at various points about 
the structure. The response data gathered is then analyzed by 
making a transformation from physical (spacial) coordinates to 
modal coordinates and the frequency response characteristics 
of the system are generated mathematically by observing the 
relationships between the forcing functions, response data, 
and known modal parameters. This allows the generation of 
mode shapes at resonant frequencies of the structure that can 
be compared to the frequencies and mode shapes provided by the 
analytical model. 

The mode shapes generated from the physical test can be 
mathematically compared to those from analysis by a number of 
methods, with one of the more popular being the cross 
orthogonality relation 

[c]^»[ W r W d- 1 ) 

where [C] Aj is the cross orthogonality matrix correlating the 
i-th analytical mode shape with the j-th test mode shape [2]. 
Matrices [<#> x ] and [0 T ] represent the assemblage of analytical 
and test acquired mode shapes, respectively, while [M A ] is the 
analytically derived mass matrix. 



3 


Under ideal correlation where the test and analysis mode 
shapes correlate identically, the resulting cross 
orthogonality matrix appears as a unit diagonal matrix. Under 
more realistic situations where the degree of correlation is 
less than ideal , the degree of variation between the two sets 
of mode shapes is identified by the density and magnitude of 
the off-diagonal terms appearing in the cross orthogonality 
matrix. Generally, satisfactory correlation is obtained if 
the cross orthogonality matrix possesses no off-diagonal terms 
exceeding 0.20 indicating twenty percent miscorrelation 
between the i-th analytical and j-th test mode shapes. 

Unfortunately, the process of modal testing can be quite 
expensive and time consuming. It is not uncommon for the 
modal test to consume many hundreds or thousands of man-hours 
to set up, conduct, and post-process the data. It also 
requires the use of specialized digital processing equipment 
and software that may not always be available to the engineer. 
This can make the correlation process somewhat prohibitive for 
smaller scale engineering projects. 

The process of correlating a finite element model can be 
optimized if it may be done without the testing required with 
the modal technique. In the case of space flight hardware, 
this optimization could be realized if the frequency response 
information needed from the test article could be provided by 
means of the random vibration data that is normally acquired 
during the flight qualification process. Random vibration 
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testing is usually required in the engineering process to 
prove the system's structural and functional survivability 
when subjected to the anticipated flight environment. This is 
in sharp contrast to the modal technique which is necessitated 
by the desire to validate the analytical assumptions. 
Consequently if the engineering assumptions can be quantified 
by the random vibration test, additional resources need not be 
allocated for a separate modal survey. 

Normally, the random vibration response data is not 
accepted as valid correlation input since it only shows the 
existence of resonance but in its raw form does not provide 
any insight into the mode shape present at the resonant 
frequency. In essence, the resonant frequency indicated by 
the random response data may ideally match a frequency 
provided by analysis yet contain a completely different mode 
of vibration. 

1 . 2 Objective 

Consider a scenario in which the random vibration data 
gathered at a discrete number of spatial locations on the test 
article is reduced to provide peak response and amplification 
data for a limited number of primary frequencies indicated by 
the test data and assumed to correspond with predetermined 
analytical frequencies. These analytical frequencies at which 
correlation is sought could be selected from the free 
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vibration analysis to represent the fundamental system 
frequencies . This can easily be accomplished by selecting 
modes in which there is significant mass participation. 
Assuming that for any of these predetermined frequencies if 
the structure were subjected to static loads equivalent to 
those generated by the flight random vibration environment, 
some sort of deflected shape could be expected. It stands to 
reason therefore that since the random vibration test provides 
local accelerations at discrete spacial locations and 
frequencies the analyst should be able to surmise a deflected 
shape for any of the resonant frequencies chosen from the test 
data. 

Given Power Spectral Density (PSD) output from a random 
vibration test, peak response and amplification values can be 
extracted for any number of selected frequencies. By use of 
basic dynamic relationships, this response and amplification 
data can be used to generate equivalent local accelerations 
corresponding to the accelerometer locations and directions 
monitored during the testing event. Assuming that these 
equivalent local accelerations possess the proper orientations 
and magnitudes, when applied to the FEM and a static analysis 
performed, a static deformation is produced resembling the 
mode shape expected at the frequency from which the response 
values were taken. This static deformation can then be mass 
normalized in the same manner as an eigenvector would be to 
form a static or 'pseudo' mode shape. This static mode shape 
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can then be used in the cross orthogonality correlation, 
taking the place of the test acquired mode shape that would 
have been generated via the modal technique. By carrying out 
this procedure, a set of static mode shapes can be generated 
from the frequency response characteristics of the physical 
model without the need for any additional testing above that 
done in the normal flight qualification process. 

JU3 Scope 

This thesis provides a detailed method for accomplishing 
the dynamic to static transformation which allows for the 
generation of the static or 'pseudo' mode shapes. 

The equations governing the correlation process, from the 
basic eigenproblem through the cross orthogonality check are 
outlined to examine their relevance to the process. A 
discussion is also presented on the preparation of both the 
analytical model and random vibration test data. This 
includes the reduction of the analytical model to those same 
degrees of freedom observed in the vibration test as well as 
the reduction of the random vibration data to the parameters 
required in the calculations. The calculation of the static 
mode shapes, their mass normalization, and comparison with the 
mode shapes provided by the free vibration analysis is 
presented in detail. 

Lastly, the correlation process is demonstrated in full 
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for both a simplistic cantilevered beam model as well as for 
an actual article of space flight hardware. The results of 
both examples indicate that the process generates accurate 
system mode shapes and provides a greatly abbreviated 
correlation method. 


CHAPTER II 


THEORETICAL DEVELOPMENT 
2_JL Fundamental Dynamic Relationships 

The basic response characteristics (for both free and 
forced vibration) of a mechanical system can be expressed by 
means of the frequency response function (FRF). This 
mathematical relation contains the mass, stiffness, and 
damping properties of the system and can be used to evaluate 
the frequency of vibration as well as the relative 
displacements of the points within the system. In terms of 
the spacial model, the mechanical system is described by mass, 
stiffness, and damping properties 

[«]ix}+[B]ba+[if]i?d=iF(t)} (2.i) 

Considering the response to be a time varied function and 
with the presence of damping the complete solution to be 
complex, the FRF of the mechanical system can be written in 
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terms of the modal properties of frequencies and eigenvectors 
(mode shapes). 

x { t) =xe iat (2.2) 


( -to 2 [M] +io> [B] + [K] ) xe iwt =Fe iwt =F( t) ( 2 • 3 ) 


a^vto) 



( r^> 7 ) [Aj 
(F r -t0 2 Af r ) +i (<0B r ) 


(2.4) 


The preceding FRF represents a single term in the 
complete FRF matrix that describes all modal properties of the 
system. Ewins [1] provides a detailed development of the FRF 
and its relevance in vibratibn testing. 

The response of the system to random excitation is a 
function of the total FRF and the source excitation. If this 
input excitation is a time varied function of acceleration and 
is considered to be a random process with a uniform spectral 
density, S x ( o>) , then the spectral density response, S 0 (<•>), may 
be represented by the following function 

S 0 (w) = |a («) | 2 S jr (w) (2.5) 

Given the Power Spectral Density (PSD) function above 
[ 3 ], it is evident that the response of a mechanical system to 
a random excitation is dependent upon the full complement of 
frequency response functions used to describe the mechanical 
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system. Realization of this dependency enables the 
transformation from random response characteristics back to 
static response characteristics as proposed in the correlation 
process. In the modal method, the frequency response function 
is operated on directly to provide the system mode shapes. 
The process to be developed uses the response as provided by 
Equation 2.5 to generate the system mode shapes. It is 
important to note, however, that the response provided by the 
random response function contains all of those same frequency 
response characteristics present in the modal method. 

Although the origin, development, and theories behind the 
frequency response and power spectral density functions are 
very important in the analysis of dynamic systems, their 
direct applications have very little impact on the development 
of the process at hand. It is left up to the reader to become 
familiar and comfortable with the mathematical origins of 
these topics. 


2.2 Preparation of Model 

The correlation process deals with three varying levels 
of fidelity representing the same mechanical system. The 
physical test specimen represents a continuum system while the 
analytical model represents the continuum system through 
significantly reduced degrees of freedom. Thirdly, the 
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testing process gathers response data at a discrete number of 
points (and directions) that represents still yet a further 
reduction in the active degrees of freedom when compared with 
the already reduced analytical model. Consequently in working 
with the analytical and test models there needs to exist a 
correspondence for those active degrees of freedom coincident 
to both. Generally this means that the positioning of 
accelerometers on the physical specimen are judiciously 
selected so that the response measurements can accurately 
represent the dominant response characteristics of the 
structure. The positioning of these accelerometers can be 
greatly aided by studying the predictions provided by the 
analytical model. The analysis can provide great insight into 
the system modes containing significant amounts of mass 
participation as well as indicating spacial locations where 
maximum dynamic deflections occur. Using this information to 
place the accelerometers increases the likelihood of observing 
the peak response characteristics during the testing event. 

However, in the reverse process of applying measured data 
to the analytical model in performing the correlation process , 
it is also desirable to reduce the many degrees of freedom 
present in the model to those same degrees of freedom which 
were measured in test. Doing so allows a one-to-one 
correlation between the analytical and physical properties. 
When using MSC/NASTRAN for the analysis, this reduction in 
active degrees of freedom can be accomplished through Guyan 
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(static condensation), Generalized Dynamic Reduction, or 
component modal synthesis. The Guyan reduction method will be 
used throughout this work. In the terminology of MSC/NASTRAN 
[4] the degrees of freedom in the displacement set (x f ) are 
partitioned into the active set (x.) and the omitted set {x 0 }. 
The active set is that set remaining after condensation on 
which the analysis is performed. Partitioning of the spacial 
system is performed such that properties of the full structure 
are maintained after condensation. In short, the reduction of 
the free vibration problem is done in the following manner 

[AT ff ] lx}+ [£ ff ] txk [K ff ] iXfMO) (2.5) 

K ~ K& ° W - j°l (2.6) 

Koa *«JW W 


[K aa ]\xJ+[K ao ]bcJ={0\ (2.7) 

[K^bcJ+lK'JbcJ^O} ( 2 . 8 ) 

[G oa ]=-[K 00 ]-i[K oa ] (2.9) 

ix^lG^ixJ ( 2 . 10 ) 
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( 2 . 11 ) 


[I] ' 

[Goa] . 


T 


[ -CD 2 [M ff ] +1(0 [B f f] + [K ff ] ] 


' [I] 

t^oa] 


Oc^Mo) 


( 2 . 12 ) 


which results in the final set of equations for the reduced 
model (system) 

[ -o> 2 [M aa ] +io> [B aa ] + [K aa ] } ixj={0\ (2.13) 

For a forced vibration problem, the static condensation 
presented above takes some liberties with the time dependency 
of the force vectors involved. But for the free vibration 
problems presented in the correlation problem, these 
assumptions pose no degradation in the solution process. 
Therefore, the Guyan Reduction methods used by MSC/NASTRAN can 
be used without concern. Alternately, Generalized Dynamic 
Reduction may also be employed (which retains all factors in 
the force vectors) with the same results being obtained. 
Regardless of the method employed, either method allows for 
the reduction from the full set of active degrees of freedom 
to a set corresponding with the test measurement points and 

directions . 

During the testing event , acceleration response data are 
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gathered and processed to provide a set of PSD curves for each 
accelerometer on the test specimen as well as those 
accelerometers mounted on fixturing which provide a control 
feedback loop to the excitation equipment. These PSD curves 
generally appear as shown in Figure 2.1, one existing for each 
DOF recorded, providing acceleration response data in terms of 
acceleration squared per unit frequency ( g 2 /Hz ) versus cyclic 
frequency (Hz). When biaxial or triaxial accelerometers are 
used, two or three such PSD curves appear for the same 
discrete spacial location at which the accelerometer is 
attached. This enables the engineer to observe response in up 
to three translational directions at any given point on the 
structure. In general, only translational accelerations are 
measured since the measurement, processing, and application of 
rotational accelerations is seldom used. 

The random response data can provide a means of 
identifying principle modes of vibration. Each set of data 
(in reference to a given excitation) can be evaluated to 
indicate principle modes excited by the input energy. 
Generally, the response data can indicate a primary mode of 
vibration by numerous appearances of distinct peaks at the 
same frequency. Once a number of these modes have been 
identified as being primary system modes or local modes of 
interest, the peak values for each can be recorded from the 
set of PSD curves. Amplification data can be calculated from 
the response data at each of these frequencies by dividing the 
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FIGURE 2.1 

Typical Response PSD Plot 
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response PSD by the input PSD (from the control 
accelerometers). For a given DOF and frequency the associated 
amplification, or quality factor, may be calculated as 


or 


PSD reS p 
PSL ^ input 


(2.14) 


2.3 Calculation of Static M ode Shapes 

Once the test data has been reduced to a set of 
frequencies, peak PSD and amplification values, the 
corresponding static or r pseudo ' mode shapes can be calculated 
using the analytical model. 

With the reduction of the analytical model to the degrees 
of freedom measured in test, as demonstrated earlier, the 
apparent mass and stiffness of the FEM are given by matrices 
[M„] and [ K, a ] * In order to produce statically deformed 
shapes from this mass and stiffness, the peak PSD and 
amplification values must be transformed to equivalent local 
accelerations so that the static system of equations may be 
solved 

[JC„] W,- [M.JIa), (2.15) 

where {xh represents the static deflection, or static mode 
shape, resulting from the application of the equivalent local 
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accelerations (a^. The vector of equivalent local 
accelerations (aK can be assembled from the individual peak 
PSD and amplification values for each DOF present. The 
transformation from dynamic response to equivalent local 
acceleration is performed by use of Mile's equation [5] 

a i =±3.0g^(^)Q 1 PSD 1 f (2.16) 

where aj = equivalent local acceleration 

g = acceleration due to gravity 
Q a = amplification factor from test 
PSD* = peak PSD value from test 
f = frequency (cyclic) of interest 
i = degree of freedom 

This local acceleration term, a lf represents the 
equivalent acceleration obtained for a single degree of 
freedom on the test article at a single frequency point. If 
the article is instrumented with five triaxial accelerometers, 
then there will be fifteen such local accelerations calculated 
and assembled to form a single acceleration vector, {ajj. 
Likewise, if four frequencies are chosen for correlation, 
there will be four acceleration vectors, {a} 1# assembled from 
the test data. 

The PSD and amplification factors involved come from test 
data at a specific frequency value that may or may not be 
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identical to the frequency indicated by analysis but at which 
the correlation is sought. To prevent introduction of error 
from the test sequence (in reference to the frequency 
discrepancies) it is preferable to use the frequency (f) as 
provided by analysis. In other words, the correlation is 
always going to be less than ideal since the test and analysis 
models will never be identical. The intended purpose of the 
correlation process is to compare the mode shapes of the two 
and therefore it is desirable to avoid introducing additional 
error into the calculations. A comparison of test and 
analysis frequencies can be made on a much smaller scale as 
will be demonstrated later. For this reason, the analytical 
frequencies are used in the calculation of the eq uivalent 
local acc elerations. 

One last step is required in the preparation of the test 
data. The response data coming from test is assumed to be in 
the proper relative magnitudes but contains no information 
concerning the orientation, or direction, of the acceleration. 
Some insight is needed to orient the applied accelerations so 
that they are meaningful in the solution of the static system 
of equations. The orientation of each acceleration can be 
taken from the corresponding degree of freedom of the 
eigenvector as predicted by the normal modes analysis. The 
analytical eigenvector for each mode to be correlated can 
easily be reduced to a directional vector containing positive 
or negative values of unity to give the accompanying 
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acceleration vector the proper orientation at the reduced 
nodal degrees of freedom. 

From the static equations, the deformed shape, {x} A 
resulting from the application of the equivalent local 
accelerations is performed by pre-multiplying both sides of 
the static equation by the inverse of the reduced stiffness 
matrix [K aa ]. 

[K„] ■’ [K.J to,- [K„] - 1 [M„] (a!, (2.17) 

to,-[J CJ " l [M„] to, (2.18) 


This statically deformed shape is not unlike a raw 
eigenvector having unbounded magnitudes of its components , and 
may appear somewhat meaningless when first observed. However, 
just as in an eigenanalysis, this deformed shape can be mass 
normalized to provide the static or 9 pseudo' mode shape. 





(2.19) 


Ideally, it would appear that the series of equations 
shown above could be applied consecutively for all the vectors 
of equivalent local accelerations representing the modes 
desired to be correlated. But in fact if this were carried 
out consecutively for a number of modes, it becomes evident 
that the lower order modes have a large influence on the 
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subsequent inodes derived. The higher order inodes exhibit 
strong content of the lower inodes, with the lowest inode 
exhibiting the greatest influence. 

To carry out the correlation process for inodes other than 
the lowest analytical mode presented, the influence of each 
successive mode must be removed from the active set of 
equations. This is easily performed by applying classical 
matrix deflation techniques [6] to the system of equations. 
Letting 

[23] i =[Ji - ]2 1 [M a .] i (2-20) 


and 


[L>] 


( 2 . 21 ) 


Then by classical matrix deflation techniques, the 
apparent mass and stiffness of the system excluding the 
effects of vectors just iterated upon are established by 


tf>J i+a =[D]i--^( 




( 2 . 22 ) 


where \ and {<f>} x are the eigenvalue and eigenvector as 
provided by analysis, not test. Using the analytical values 
again keeps the active mass and stiffness from becoming 
contaminated with the error introduced by test. Also note 
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that to prevent the introduction of error, the eigenvalue and 
eigenvector taken from analysis must be accurate enough to 
prevent affecting downstream calculations. 

If the analytical eigenvector is mass normalized as part 
of its original processing, the denominator in the right side 
of the equation takes on the value of unity by definition. 
This allows the further simplification to 

[D] [D] r ^}#I[ M aa ] (2.23) 


The solution of the second deflected shape in the series is 
then 

W it r[o] i+ M^ ( 2 . 24 ) 

This deformed shape is also mass normalized with respect 
to the analytical reduced mass matrix to complete the 
iteration and all subsequent static mode shapes are iterated 
upon in the same manner . 

It is important to reiterate that each successive mode 
shape appearing in the analytical results must be swept from 
the mass and stiffness matrices for the calculation of higher 
order modes, even if that particular mode is not one to be 
correlated. That is to say that if it is sought to correlate 
the first, third, and fourth analysis modes, the second mode 
must be swept from the system before the third and fourth 
eigenvectors are calculated. Otherwise the presence of the 
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second mode will heavily influence the latter static mode 
shapes . 


ZjlA Correlating the Mode Shapes 

When all of the static mode shapes desired for 
correlation have been calculated, the degree at which they 
match the analytical mode shapes can be evaluated using any 
number of methods, e.g. Modal Assurance Criterion, Cross 
Orthogonality Criterion, etc. For the sake of demonstration 
purposes, the cross orthogonality method has been chosen in 
this thesis work. The cross orthogonality method can be 
expressed by 

[C] =[<!>*] r [Af a J [* T ] (2.25) 

Where [0 A ] is a matrix formed by assembling the eigenvectors 
from analysis and [0 T ] is the matrix formed by assembling the 
static or ' pseudo ' mode shapes as calculated by the process 
outlined above. The cross orthogonality matrix is a square 
matrix where the diagonal terms represent the orthogonality 
condition between each of the analytical and test 
eigenvectors. If the i-th analysis eigenvector is identical 
to the i-th test (pseudo) eigenvector, the i-th diagonal term 
in [C] would be unity. 
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= (2.26) 

Therefore anything less than unity on the diagonal indicates 
some level of skew between the analysis and test eigenvectors. 

The off-diagonal terms represent the orthogonality 
condition between the i-th analytical and j'-th test 
eigenvectors. By the principal of orthogonality, two 
eigenvectors from the same eigenproblem must be perpendicular, 
or orthogonal, to each other. Mathematically 

(2 - 27) 

Therefore for ideal correlation between the analytical 
and test models, the cross orthogonality matrix would appear 
as a unit diagonal matrix. Realistically, where the 
correlation between the analytical and test models is not 
ideal the diagonal terms will become less than unity while the 
off-diagonal terms will increase from zero. If the 

correlation experienced is less than acceptable, the finite 
element model can be 'tuned' to gain more accurate 
distributions of mass and stiffness. Iterations on the FEM 
can be followed by recalculation of the static modes shapes 
and the cross orthogonality matrix. The vibration testing 
need only be conducted once since the same data is used for 
each iteration on the finite element model. Usually, 
satisfactory correlation has been achieved when the absolute 
values of diagonal terms are greater than 0.90 and the 



24 


absolute values of the off-diagonal terms are limited to 0.20 
or less. These limits are usually applied for those modes 
deemed principal system modes or for local modes of particular 
interest. 


2.4 Automating the Process 

A Direct Matrix Abstraction Program (DMAP) alter has been 
written for Version 67.5 of MSC/NASTRAN [7] to carry out the 
sequence of operations presented above. This DMAP alter uses 
the analytical model, once reduced with Guyan or Generalized 
Dynamic Reduction, to calculate the static mode shapes in the 
same run sequence as a free vibration analysis is performed. 
The DMAP statements are compiled in the MODERS module where 
MSC/NASTRAN performs the real eigenvalue analysis. The user 
is required to input matrices assembled from the vectors of 
peak PSD and amplification values (Q) as well as a listing of 
the analytical modes to be correlated. Peak PSD and 
amplification values are entered via the DMI bulk data entry 
cards (PSD and QUAL) while the list of modes to be correlated 
are entered via a DTI entry (MODELST). The user is also 
required to input several constant parameters to describe the 
acceleration due to gravity (GRAV) , maximum number of modes 
predicted by analysis (NMAX) , number of modes to be matched 
(NMODES) , and number of reduced degrees of freedom (NTEST) . 
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The latter three values are used for partitioning the various 
matrices assembled internal to MSC/NASTRAN. The user is also 
required to disable the internal resequencing routine within 
MSC/NASTRAN so that no transformation between the internal and 
external node numbering schemes is implemented. This helps to 
prevent application of the external accelerations to wrong 
internal degrees of freedom. 

The DMAP alter recovers the reduced analytical mass and 
stiffness matrices, M aa and K aa , as well as the eigenvalue and 
eigenvector matrices from the modal analysis. These 
eigenvalue and eigenvector matrices are partitioned and 
reduced for the modes to be correlated as indicated in the DTI 
MODELST entry. The reduced eigenvector matrix is partitioned 
term by term and the directional orientation of each DOF is 
determined. A matrix (NDIR) is then assembled using these 
directional orientations. The test data is then read and the 
equivalent local acceleration matrix is formed by using Mile's 
relation with the information provided. 

The static mode shapes are calculated and mass normalized 
as demonstrated earlier. After each static mode shape is 
calculated, the influence of the corresponding analytical mode 
shape is swept from the reduced mass and stiffness matrices. 

intermediate modes not chosen for correlation are also 
swept from the system before subsequent higher order static 
mode shapes are calculated. The resulting matrix of static 
mode shapes (PHITST) is automatically printed for the user. 
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Lastly, the reduced analytical eigenvector, reduced analytical 
mass, and static mode shape matrices are used to perform the 
system cross orthogonality calculation. The resulting cross 
orthogonality matrix (COC) is also printed for the user. 

A complete listing of the DMAP alter is given in Appendix 
A. 



CHAPTER III 


CORRELATION PROCESS 


3, 1 Demonstration I 


For purposes of developing the correlation process, a 
simple cantilevered beam model was used as the baseline 
system. The analytical model is a cantilevered beam as shown 
in Figure 3.1 consisting of four (4) nodes with two degrees of 
freedom per node, one translational and one rotational. The 
nodes are connected by three Bernoulli-Euler beam elements 
with consistent mass and stiffness [8] as given below 
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E = 70.0 MPa 
V = 0.33 
p = 2700 kg/m 3 

I = 5.0 X 10 -10 m 4 
A = 6.0 X IQ" 5 m 2 


FIGURE 3.1 

Cantilevered Beam Analytical Model 
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12 6 L -12 6 L 

r ri./M) 6L 4i2 - 6i 2L2 (3.2) 

2 L 3 -12 -6L 12 -6L 

6L 2L 2 -6Z, 4L 2 


Substituting the geometric properties of the beam, 
assembling the connectivity, and condensing out degrees of 
freedom associated with the cantilevered restraint, the mass 
and stiffness matrices appear as follows 

312.00 0.00 54.00 -1.30 0.00 0.00 

0.00 0.08 1.30 -0.03 0.00 0.00 

54.00 1.30 312.00 0.00 54.00 -1.30 

-1.30 -0.03 0.00 0.08 1.30 -0.03 

0.00 0.00 54.00 1.30 156.00 -2.20 

0.00 0.00 -1.30 -0.03 -2.20 0.04 

(3.3) 

840.00 0.00 -420.00 21.00 0.00 0.00 

0.00 2.80 -21.00 0.70 0.00 0.00 

-420.00 -21.00 840.00 0.00 -420.00 21.00 

21.00 0.70 0.00 2.80 -21.00 0.70 

0.00 0.00 -420.00 -21.00 420.00 -21.00 

0.00 0.00 21.00 0.70 -21.00 1.40 

(3.4) 

For f ree-vibration , the basic eigenproblem can be given as 

[M aa ] txk [K aa ] ixMO} (3.5) 



[M aa ] =3 . 8 5x1 O' 5 
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Assuming the solution is harmonic in nature, 

x i =Xcos (<ot-a) (3.6) 

i i = -ti)Xsin(Q t-a) (3.7) 

£ i = -u> 2 Xcos {ot-a) (3.8) 

Therefore the complete solution is 

-u 2 Xcos (wt-a) [M aa \ +^cos(«t-a) [K aa ] =(0f (3.9) 

[[K aa ] -w 2 [M aa ] jOdcos (ot-a) =(o} (3.10) 

[[K aa ]-<* 2 [M aa ]]M={ 0} (3.11) 

For which the nontrivial solution is 

det[ [K aa ] - 0 > 2 [AfaJ ]=0 (3.12) 


Using inverse iteration [6] as the solution scheme and the 
Rayleigh quotient to force convergence of the eigenvalue 

[D] x = [ K aa ] i 1 [M aa ] 1 


( 3 . 13 ) 
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<v s *i }= u>i iW 


(3.14) 


k s+1 }= 




v /tv s+1 F[Af aa ](v s+1 ) 


(3.15) 


and the associated eigenvalue is 


fi+1 (v s+1 ) [M aa ] 


(3.16) 


Upon satisfactory convergence of the eigenvalue, the 
final eigenvector can be normalized with respect to the system 
mass (mass normalization). 



^v s+1 F[Af aa ] (v s+1 ) 


(3.17) 


Subsequent eigenvectors and eigenvalues are found in the 
same manner once the system is shifted to remove the influence 
of the present mode. Therefore 

[D] i+1 = [D] [M^] (3.18) 

Performing the series of iterations on the analytical 
model produces the following results 



32 




' 329.8 

1 . 304e4 1 . 041e5 

5 . 278e5 

1 . 87 0 e6 

7.430 e6 



2.89 

18.36 

51.35 

115.6 

217 .6 

433 . 8 

■ V 
€ 


1.5021 

-5.3845 

6 .7827 

2.4709 

2.7984 

2.3954 


27 .358 

-53.667 -50.501 

-315.713 

521.238 

245.673 

£i 


4.9630 

-3.8654 -5.9809 

1.1416 

-4.9303 

3.7505 


39.591 

90.048 -38.548 

323.470 

342.850 

690. 579 



9 . 0740 

9.1282 

9.1159 

-9.8090 

-10.6002 

19.9162 



41.635 

145.594 241.534 

-418.042 

-702.697 

1990.912 


(3.19) 


The validity of the eigenvectors calculated can be 
reinforced by examining the orthogonality conditions present. 
By definition of mass orthonormality, the eigenvectors 
associated with distinct eigenvalues must satisfy the 
following relation 


[<M r [w] W = in 


(3.20) 


where [I] is the real identity matrix. Carrying out the 
operation 


[<w T [M aa ) [<w 4 


1.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

0.0000 


0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 


0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 


0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 


(3.21) 


10- 7 


All errors in the orthogonality test were on the order of 
indicating that the eigenvectors satisfy the 
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orthonormality constraint. 

The test or physical system for this demonstration 
problem was simulated with a higher fidelity model (to imitate 
a continuum system) . An MSC/NASTRAN model was created with 
sixty-one nodes, again with two degrees of freedom per node, 
and sixty CBAR structural elements. Coupled mass formulation 
was employed as in the analytical model. This model is shown 
in Figure 3.2. 

Initially, a normal modes analysis was performed on this 
'test' model to compare the resulting frequencies with those 
from the lower fidelity analytical model. The comparison of 
eigenvalues and cyclic frequencies is shown in Table I. 

As can be seen in the comparison, the lower fidelity 
analysis model is significantly over-stiffened in the higher 
modes. The stiffening effect of the analytical model is shown 
graphically in Figure 3.3. The systematic deviation from the 
ideal line shows the uniform mathematical stiffening created 
by the dramatic reduction in degrees of freedom associated 
with the analytical model. 

The eigenanalysis of the 'test' model was followed with 
a random vibration analysis. The original model was augmented 
by attaching a large (5000 kg) mass to the support node 
applying a random acceleration excitation (Figure 3.4) to this 
mass. The excitation was applied in the translational (Y) 
direction to force bending response in the beam. 
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with Input PSD Applied 


FIGURE 3.2 

High Fidelity Cantilevered Beam Model 
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TABLE I. 

Comparison of Analytical and Test Eigenproblem Results 


Mode 

Eigenvalue 

Frequency (Hz) 

Anal. 

Test 

Anal. 

Test 

%Diff 

1 

329.8 

329.7 

2.89 

2.89 

0.00 

2 

13305.3 

12950.1 

18.36 

18.11 

-1.42 

3 

104. le3 

101. 5e3 

51.35 

50.71 

-1.25 

4 

527. 8e3 

389. 9e3 

115.6 

99.37 

-14.0 

5 

186. 9e4 

106. 5e4 

217.6 

164.3 

-24.5 

6 

743. 0e4 

237. 8e4 

433.8 

245.4 

-43.4 




































Test Frequency, Hz 
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Analysis Frequency, Hz 


FIGURE 3.3 

Comparison of Test and Analytical Frequencies - Cantilevered 

Beam Model 
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Frequency, Hz 


FIGURE 3.4 

Excitation PSD at Large Mass 
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Structural damping equal to ten percent of critical 
damping (Q = 5) was included in the structure and was held 
constant throughout the frequency range of zero to two 
thousand Hertz. The translational response Power Spectral 
Densities were recovered at the three nodes corresponding to 
the node locations present in the analysis model (excluding 
the support node). These PSD plots are shown in Figures 3.5, 
3 . 6 and 3.7. The peak response values at resonant frequencies 
were recorded from the analysis and are summarized in Table 

By dividing each of the response predictions by the input 
(control) acceleration, the response amplification values can 
likewise be tabulated for each active degree of freedom. This 
set of amplification factors is shown in Table III. 

Referring to the matrix of eigenvectors from the 
analytical model (see Eq. 3.19) the directional attributes of 
each DOF in each mode can be determined. The nodal 
directional coefficients are shown in Table IV. 

Using Mile's relation as presented earlier, the 
combination of peak response, amplification, and orientation 
are used to establish the set of equivalent local 
accelerations. Note that the frequencies established by the 
analytical model are used in the calculation of these 
accelerations. These accelerations are shown in Table V. 

Applying these equivalent local accelerations to the 
analytical model using the process outlined earlier and mass 
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HIOJ FIDELITY CANTILEVERED BEAM MODEL 


FIGURE 3.7 

Response PSD at Full Span 
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TABLE II. 


Peak Response Values 


Peak PSD Values (g^Hz) 


Freq/DOF 


Input 


21T2 


21R3 


41T2 


41R3 


61T2 


61R3 


2.89 


2.890e-4 


8.689e-4 


0.0 


5.595e-3 


0.0 


1.799e-2 


0.0 


18.11 


1.881e-3 


1.029e-2 


0.0 


5.950e-3 


0.0 


3.012e-2 


0.0 


50.71 


1.300e-2 


4.675e-2 


0.0 


3.365e-2 


0.0 


8.031e-2 


0.0 


99.37 


1.500e-2 


4.908e-3 


0.0 


2.146e-3 


0.0 


4.531e-2 


0.0 


164.3 


1.288e-2 


9.372e-3 


0.0 


6.995e-3 


0.0 


2.477e-2 


0.0 


245.4 


5.290e-3 


1.919H 


0.0 


6.339e-3 


0.0 


1.243e-2 


0.0 
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TABLE III. 

Response Amplification Values 
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TABLE IV. 

Acceleration Orientations 


Freq/ 

Acceleration Orientations 

DOF 

2.89 

18.11 

50.71 

99.37 

164.3 

245.4 

21T2 

1.0 

-1.0 

1.0 

1.0 

1.0 

1.0 

21R3 

1.0 

•1.0 

-1.0 

-1.0 

1.0 

1.0 

41T2 

1.0 

-1.0 

-1.0 

1.0 

-1.0 

1.0 

41R3 

1.0 

1.0 

-1.0 

1.0 

1.0 

1.0 

61T2 

1.0 

1.0 

1.0 

-1.0 

-1.0 

1.0 

61R3 

1.0 

1.0 

1.0 

-1,0 

-1.0 

Em 
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TABLE V. 

Equivalent Local Accelerations 
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normalizing each resulting vector with respect to the reduced 
analytical mass matrix, the following static mode shapes can 
be derived. 


‘ V 

fi 


■ 329.8 
2 . 89 
1.473 
27 . 046 
4.940 
40 . 083 
9 . 128 
42.523 


1 . 331 e 4 
18.36 
- 5.322 
- 52.912 
- 3.934 
88.439 
9.251 
149 . 021 


1 . 041 e 5 
51.35 
6 . 819 
- 48.600 
- 6 . 018 
- 34.436 
8.940 
230.603 


5 . 27 8 e 5 
115.6 
2.306 
- 353.950 
1.962 
331.998 
- 7 . 07 0 
- 189 .772 


1 . 87 0 e 6 
217 .6 
3.525 
576.020 
- 3 . 009 
591.330 
- 1.709 
156 . 038 


7.430 e 6 
433 . 8 
- 2.395 
- 245.673 
- 3.751 
- 690 . 579 
- 19.916 
- 1990 . 91 . 


( 3 . 22 ) 


As demonstrated earlier, the orthogonality condition of 
the eigenvectors calculated can be evaluated using Equation 


3 . 20 . 


[<M r [M aa ] [* r ] - 


1.000 0.006 
0.006 1.000 
0.000 0.010 
0.000 - 0.004 
0.000 0.000 
0.000 0.000 


0.000 

0.000 

0.010 

- 0.004 

1.000 

0.006 

0.006 

1.000 

0.007 

- 0.074 

0.001 

- 0.073 


0.000 

0.000 

0 . 000 

0 . 000 

0.007 

0 . 001 

- 0.074 

- 0.073 

1.000 

- 0.402 

- 0.402 

1 . 000 


( 3 . 23 ) 


The orthogonality results shown above indicate that 
static eigenvectors calculated do experience some error, 
particularly with respect to the orthogonality between the 
fifth and sixth modes. The remaining modes appear to modestly 
satisfy the orthogonality constraint. The fifth and sixth 
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modes appear to be fairly orthogonal to the other four modes 
but do experience some degree of skew between each other. 

Using the set of static eigenvectors and considering them 
to be the test acquired mode shapes, the cross orthogonality 
check can be performed. 


1.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.006 

1.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.010 

1.000 

0.000 

0.000 

0.000 

0.000 

- 0 . 004 

0.007 

0.991 

0.000 

0.000 

0.000 

0 . 000 

0.008 

- 0.113 

0.916 

0.000 

0 . 000 

0 . 000 

- 0.001 

0.073 

0.402 

- 1.000 


( 3 . 24 ) 


As can be seen, five out of the six modes correlate 
nearly identically. The fifth mode correlates to a lesser 
degree with the analytical modes as exhibited by the decreased 
diagonal term and the larger off-diagonal term. This 
indicates some level of skew from the adjacent eigenvectors. 
If the generalized stiffness for this mode is calculated 

[ K a ] ={ 4 >} T [ K aa ]{ tf =\ ( 3 . 25 ) 


[ tf e ] 5 ={< M 5 T [•*■«] {< M 5 = 2 . 776 x 10 6 > 1 . 87 0 X 1 0 6 ( 3 . 26 ) 


This shows that the fifth mode is predicted to be much stiff er 
than that predicted by analysis. By comparison, if the 
generalized stiffness of the sixth mode is calculated in the 
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same manner 

[jy 6 ={<j> 3 ) 6 r [K aa ] =7 . 43 024X1 0 6 «7 . 43 02 OxlO 6 (3.27) 


This sixth mode shape is nearly identical to that predicted by 
analysis. 

Overall this model would be considered to be well 
correlated even with the presence of the larger off-diagonal 
term relating to the fifth mode. Nearly ninety-one percent of 
the total effective mass is present in the first four modes 
correlated, with only four percent occurring in the fifth and 
sixth modes combined. This indicates that the principle 
vibratory modes are well represented by the static mode shapes 
that were generated by process. 

Formatted input and output from the MSC/NASTRAN analyses 
of the high fidelity beam model are given in Appendix B. 


3.2 Demonstration II 

Although the process has been demonstrated on a small 
scale theoretical model, the entire process is intended to be 
used during the engineering development cycle. Therefore, the 
correlation process will also be applied to an actual space 
flight structure for which the finite element modeling and 
random vibration test sequences have already been performed. 
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The structure in question is a NASA space flight structure 
that flew aboard Space Shuttle Columbia during the First U.S. 
Microgravity Laboratory (USML-1) Spacelab Mission in 1992. 
The package is approximately 19 inches wide, 36 inches tall, 
26 inches deep, and weighs approximately 300 pounds. The 
primary structural component is a large aluminum plate 
machined with integral stiffeners to form a waffleplate, or 
stiffened panel. This waffleplate supports various optical 
components and the experiment test chamber. The waffleplate 
is in turn supported by a framework of large channel sections 
running front to back at all four corners. The front and rear 
faces of the package are constructed of aluminum sheet that 
form shear panels to tie the various structural components 
together. Once mounted in the Spacelab module, the entire 
structure is supported by the front shear panel and four 
brackets attached to the primary channels at either of the 
four rear corners. Several photos of the flight system are 
shown in Figures 3 . 8 and 3.9. 

An MSC /NASTRAN finite element model was constructed 
during the engineering design process to serve as both the 
dynamic and stress model for the system. Various views of 
this finite element model are shown in Figures 3.10 and 3.11. 
The model consists of plate and offset bar elements to form 
the waffleplate structure with concentrated masses used to 
represent those components attached to the waffleplate. The 
rest of the primary structure is modeled using beam elements 
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FIGURE 3.8 

STDCE Flight Experiment Structure 
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FIGURE 3.9 

STDCE Flight Experiment Structure 
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FIGURE 3.10 

STDCE Flight Experiment Structure Finite Element Model 
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FIGURE 3.11 

STDCE Flight Experiment Structure Finite Element Model 
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to represent the primary channel sections and various 
stiffeners, while plate elements to represent the front and 
rear shear panels. A separate study of the mass properties 
showed the model to conform well with the actual structure in 
both mass and center of gravity. 

During the flight qualification process the Experiment 
Package was rigidly fixtured and subjected to random vibration 
testing at the maximum expected flight levels [9]. The input 
spectra are shown in Figure 3.12. Three separate tests were 
conducted, one in each of the three primary coordinate 
directions with single axis base excitation. Response 
acceleration measurements were taken at seven discrete 
locations on the structure. Four of the accelerometers were 
triaxial accelerometers intended to measure in-axis response 
as well as response in the two cross-axis directions. One 
accelerometer was biaxial and the last two accelerometers were 
uniaxial accelerometers mounted to the front and rear shear 
panels and intended to measure the out-of-plane panel 
response. The locations of these accelerometers are shown in 
Figure 3.13. The seven accelerometers together allowed 
response measurements in sixteen degrees of freedom. Input 
excitation was measured by four triaxial accelerometers 
attached to the test fixture, and the response of these four 
accelerometers was averaged to form the control feed-back 
scheme for the test. 

To correlate this FEM, the model was first reduced to the 
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FIGURE 3.12 

STDCE Experiment Structure Random Vibration 
Test Levels 
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same sixteen degrees of freedom at which response data were 
measured during the vibration tests. This was done by 
including those degrees of freedom corresponding to the 
accelerometers (response only) in ASET1 bulk data entries and 
employing Guyan Reduction. 

A preliminary modal analysis was conducted on the model 
to identify the principal modes of vibration. A separate 
MSC/NASTRAN alter [10] was employed in which the modal 
effective mass of each mode was calculated. This allows the 
analyst to determine which modes may be considered primary 
system modes and exclude from the correlation process local 
modes where negligible system mass participates. These 
neglected modes are often local plate modes or coupled modes 
where several modes of vibration are closely spaced. A 
frequency span from zero to two— hundred Hertz was swept for 
all modes within the range. A total of eight modes appeared 
in this range, varying from 66.5 to 186.7 Hertz. Mode shape 
plots for these modes are shown in Figures 3.14 through 3.21. 
These eight modes represent about forty-five, twenty-five, and 
forty-five percent of the total effective mass participation 
in the coordinate X, Y, and Z directions, respectively. 
Although these modes represent only a minority of the total 
effective mass, this frequency range in which they occur is 
suitable for correlating the model. Modes 1, 3, 4, 5, and 7 
have been selected for correlation based on the large 
effective mass content in each of the modes. The second and 
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FIGURE 3.14 

STDCE Experiment Package Mode 1 - 66.5 Hz 
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FIGURE 3.15 

STDCE Experiment Package Mode 2 - 71.1 Hz 
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FIGURE 3.16 

STDCE Experiment Package Mode 3 - 88.1 Hz 
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FIGURE 3.17 

STDCE Experiment Package Mode 4 - 108.3 Hz 
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FIGURE 3.18 

STDCE Experiment Package Mode 5 - 112.7 Hz 
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FIGURE 3.19 

STDCE Experiment Package Mode 6 - 122.2 Hz 
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FIGURE 3.20 

STDCE Experiment Package Mode 7 - 163.2 Hz 


65 



FIGURE 3.21 

STDCE Experiment Package Mode 8 - 186.6 


66 


sixth inodes were excluded from the correlation since they 
contain less than one percent of the total effective mass. 
The MSC/NASTRAN output listing for this effective mass 
analysis is given in Appendix C. 

To proceed with the correlation process, the modes 
determined to be primary modes now need to be identified and 
separated in the random vibration test data. This can most 
easily be accomplished by examining and sorting the test data 
based on peak PSD values and the number of degrees of freedom 
where these values appear at the same frequency. Quite often 
one particular set of data, with respect to the excitation 
axis, provides a more discernable appearance of a fundamental 
mode more than the other sets of data. The set of data from 
which the data for the mode is taken doesn't necessarily 
coincide with the primary response direction (or effective 
mass direction) of the mode since the response mode may most 
easily be excited by cross-axis input (such as with a 
torsional mode). However, these strong system modes usually 
stand out markedly over the less effective or coupled system 
modes . 

Figure 3.22 shows the averaged control and three response 
PSD plots corresponding to a single triaxial accelerometer 
affixed to the waffleplate. These plots were taken from the 
seventeen plots that make up the data set from the X— axis 
test. Here it is easy to see the existence of a resonance at 
105.5 Hertz. The peak values at this frequency have been 
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1000 20 100 1000 
Figure 3.22 

Sample Control and Response Measurements 
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indicated to illustrate how a data set is assembled. For the 
STDCE Experiment Package, the peak PSD response data for the 
five modes of interest has been reduced from the full 
complement of test data and is shown in Table VI. 

Amplification values (Table VII) are computed in the same 
manner as demonstrated earlier using the control and response 
PSD values. 

A graphical comparison of the selected test frequencies 
with the analytical frequencies is given in Figure 3.23. This 
shows the third mode (88.1 Hz) to be somewhat stiffened in the 
actual structure from what is predicted by analysis. The 
remaining modes are overly stiffened in the model as might be 
expected. These frequency differences may be the direct 
result from the reduction xn active degrees of freedom in the 
model but may also be the result of effects not accounted for 
in the modeling and linear analysis, e.g. joint friction and 
other various nonlinear effects. 

Before the PSD and amplification parameters can be used 
in the correlation process, they must be rearranged in the 
proper sequence to match the nodal numbering in the analytical 
model . This ensures that the resulting local accelerations 
are applied to the proper node and in the proper translational 
direction. With the MSC/NASTRAN nodal renumbering scheme 
deactivated, the system matrices are assembled and processed 
using the external node sequencing, relating the degrees of 
freedom with the nodes in ascending order. Therefore when the 
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TABLE VI. 


Peak Response Data at Selected Frequencies. 


ANAL. FREQ. 


ANAL. ORDER 
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TABLE VII 


Response Amplification Values at Selected Frequencies 


ANAL. FREQ. 


ANAL. ORDER 























































































fig: 
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PSD and amplification values are prepared in the DMI entries, 
the values are entered according to the corresponding node 
numbers in ascending order. Table VIII shows the node and 
accelerometer association for the Experiment Package. 

The results of the correlation run can be summarized 
in matrix form, one matrix each for the analytical mode 
shapes, calculated static mode shapes, and cross orthogonality 
results . 

-2 . 039e-2 -1.393e-2 6.638e-2 7.351e-2 3.331e-l' 

-2.974e-3 1 . 002e-l -2.645e-3 2.342e-2 -1.359e-l 

6.034e-3 -5 . 098e-2 7.210e-3 -l.024e-l 5.731e-l 

-2 . 6 83 e-2 -5.637e-2 -2.737e-2 -7.829e-2 -7.658e-2 

5.758e-l -2.385e-l 6.301e-2 -1.077e-2 -1.567e-2 

5.161e-2 1 . 320e-l -1.262e-l -2.654e-2 1.524e-l 

-1 . 416e-l -4.459 e-1 4.627e-l 2.350e-2 -1.106e-l 

-2.892e-3 1 . 215e-l 1.237e-l 6.252e-l -1.946e0 

2 . 595e-l 1.415e0 2.815e0 7.105e-l 1.224e0 

-5 . 351e-l -1 . 364ed 9.458e-l -3.395e-2 -1 .299 e-2 

-1 . 232e-3 6.458e-2 9.748e~2 6.933e~l -1.896e0 

2.029 eO -7.294e-l 1.168e0 -3.187e-l -7.265e-l 

-5.483e-l -1 . 418e<9 9.204e-l 6.940e-2 -6.299e-3 

7 . 982e-l -1.101 eO -1.111 eO 4.579 eO 1.422 eO 

-4 . 164e-l 1 . 162e-l -6.254e-2 -2.414e-2 -7.733e-l 

-1 . 881e0 -9 .126e-l 6.913e-l 4.695e-2 3.628e-l 



(3.28) 



73 


TABLE VIII. 

Accel erometer/Node Association 


Accel. 

Node/DOF 

Accel. 

Hode/DOF 

MIX 

603/T1 

H9X 

536/T1 

H2Y 

603/T2 

H10Y 

536/T2 

H3Z 

603/T3 

H11Z 

536/T3 

H4X 

744/T1 

H12Y 

224/T2 

H5Y 

744/T2 

N13Z 

224/T3 

H6Z 

744/T3 

H14X 

149/T1 

H7Y 

159/T2 

H15Y 

149/T2 

H8Y 

12/T2 

H16Z 

149/T3 
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-1 . 147e-2 -2.501e-2 4.464e-2 -1.558e-l 3.548e-l' 

-2 . 049 e-1 1 . 026e-2 -7.058e-2 5.764e-2 -1.175e-l 

2 . 179e-l -5.108e-2 7.808e-2 -2.303e-l 4.065e-l 

1 . 259e-l -6.862e-2 -1.850e-2 -1.566e-2 -6.496e-2 
-1 .116 e-1 -2.455e-l 5.173e-2 -1.052e-3 1.244e-2 

8 . 013e-2 1 . 543e-l -8.025e-2 -5.418e-2 9.557e-2 

-1 . 591e-l -5 . 262e-l 3.218e-l 2.563e-2 -2.090e-2 

-2.815e-l 6.603e-2 7.671e-2 1.157e0 -1.913e0 

3 . 011e-l 5.862e-l 2.757e0 3.358e-l 1.220 eO 

-3 . 267 e-1 -1.578 eO 9.239e-l 4.217e-2 -3.074e-2 

-1.726e-l 1.607 e-2 9.571e-2 1.205ed -1.964e0 

1.975e0 -1.029 eO 1.157 eO -2.162e-l -6.645e-l 

-1 . 826e-l -1 . 636e0 9.811e-l l.260e-l -I.108e-1 

1 .115eO -6 . 563e-l -1.452e0 3.874e0 1.498e0 

-4.581e-l 1.318e-l -7.347e-2 4.477e-4 -6.701e-l 

-1 . 426e0 -1 . O48e0 6.569e-l -2.237e-l 4.127e-l. 

(3.29) 

0.9335 0.0000 0.0000 0.0000 0.0000 

-0.1854 0.9591 0.0000 0.0000 0.0000 

0.0095 -0.2819 0.9955 0.0000 0.0000 

0.1223 0.0120 -0.0661 0.9551 0.0000 

0.1702 0.0143 -0.0045 -0.2811 0.9941 

(3.30) 


The process shows the model to correlate very well for 
the modes selected. The sparse content of off-diagonal terins 
and the relative magnitudes of the diagonal terms indicate 
acceptable correlation between the analytical model and 
physical system. The third and fifth analytical modes (second 
and fourth selected) correlate to a lesser degree than the 
other modes selected but still show satisfactory similarity. 
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Overall the analytical model may be considered to be a 
relatively good approximation of the mass and stiffness 
distributions of the continuum system. If a higher degree of 
correlation was sought, the model could be 'tuned 7 by 
adjusting mass and stiffness to provide more accurate mode 
shapes and frequencies. 



CHAPTER IV 


DISCUSSION 


4.1 Summary and Conclusions 

A method has been outlined for substantiating the free 
vibration predictions provided by finite element analysis 
using random vibration data acquired via dynamic testing. 
Power Spectral Density response data taken at a certain number 
of spacial locations on the structure and corresponding to a 
unique resonant frequency is used to generate a set of 
equivalent local accelerations that can be applied to a finite 
element model. Assuming that the finite element model is an 
accurate representation of the continuum system, and therefore 
possesses similar frequency response characteristics, these 
equivalent local accelerations deform the structure in a 
manner consistent with the mode of vibration from which they 
were derived. This statically deformed mode shape can then be 
compared with the free vibration mode shape at the 
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corresponding analytical frequency using modal techniques such 
as the cross orthogonality method. 

The process has been applied to a simple cantilevered 
beam system using two analytical models: a high fidelity 
model to simulate the continuum system and a significantly 
lower fidelity model to serve as the analytical model. 
Through random response analysis, response Power Spectral 
Density and amplification levels were predicted over a large 
frequency range to serve as data that would be acquired during 
dynamic testing. These response levels were reduced at six 
frequencies corresponding to the natural frequencies of 
vibration predicted by the lower fidelity model. The reduced 
data was used to generate sets of equivalent local 
accelerations which were applied to the lower fidelity model. 
The resulting system of static equations were solved and six 
static deformations predicted, mass normalized, and compared 
to the eigenvectors provided by the previous modal analysis. 
Five of the six static mode shapes were found to correlate to 
a very high degree with the analytical eigenvectors. 

The same correlation process was also applied to a NASA 
flight structure for which random vibration data had been 
acquired in the flight qualification process. Five primary 
system modes were chosen for correlation based on their mass 
participation content. Three of the five modes evaluated were 
found to correlate to a very high degree while the remaining 
two correlated to a lesser but still acceptable degree. The 
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correlation of these two inodes could be significantly improved 
through the process of tuning the finite element model. The 
tuning process, however, is beyond the immediate scope of this 
thesis . 

The results from both demonstrations indicate that the 
process can effectively be used to generate the mode shapes of 
the physical structure and thus correlate finite element 
results. The application of this process has the potential 
to substantially reduce the costs associated with correlating 
such models by eliminating the need for additional testing 
such as that required by classical modal techniques. All 
information required to complete the correlation process can 
be obtained during the normal flight qualification process. 


4.2 Recommendations for Further Study 

Even though the process developed here has been shown to 
provide sufficient correlation results, it has not been 
compared directly with results provided by alternate methods, 
e.g. modal testing. The full justification of merits provided 
by this method should be fully investigated by performing an 
accepted correlation procedure in parallel with the proposed 
procedure. Doing so would allow the comparison of test 
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acquired mode shapes from the two methods and also allow a 
characterization of the mathematical error that is introduced 
in the process that has been proposed. 
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$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ 

$ January 1994 

$ 

$ CROSS-ORTHOGONALITY CHECK OF ANALYTICAL MODE SHAPES 
$ USING 'PSUEDO' STATIC MODE SHAPES DERIVED FROM RANDOM 
$ VIBRATION DATA. 

$ 

$ Written by: Jeffrey A. Chambers 

$ NASA Lewis Research Center 

$ M/S 86-12 

$ 21000 Brookpark Road 

$ Cleveland, Ohio 44135 

$ (216)433-8612 

$ 

$ This DMAP alter is to be used with a Solution 103 
$ Normal Modes Analysis. It was written and tested using 
$ MSC/NASTRAN Version 67.5 for non-superelement models. 

$ 

$ Input: 

$ 

$ PHIANL - Analysis Derived Mode Shapes 

$ MANL - Analysis Derived Mass Matrix 

$ KANL - Analysis Derived Stiffness Matrix 

$ NMODES - Number of Modes to be Matched 

$ (PARAM, NMODES) 

$ NMAX - Number of Modes Provided by Analysis 

$ (PARAM, NMAX) 

$ NTEST - Number of Test Points (PARAM, NTEST) 

$ MODELST - Mode List for Correlation (DTI, MODELST) 

$ GRAV - Acceleration Due to Gravity (PARAM, GRAV) 

$ PSD - Peak Response Matrix in g A 2/Hz (DMI, PSD) 

$ QUAL - Amplification Factor Matrix (DMI, QUAL) 

$ NEWSEQ - Nodal Resequencing Deactivation 

$ (PARAM, NEWSEQ, -1) 

$ 

$ Output: 

$ 

$ COC - Cross Orthogonality Matrix 

$ PHITST - Test Derived Mode Shapes 

$ ACCEL - Nodal Acceleration Matrix 

$ NDIR - Nodal Directional Matrix 

$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ 

$ Compile MODERS processor just after OFP module. 

COMPILE SEMODES SOUIN=MSCSOU NOLIST NOREF 
ALTER 2,2$ 

TYPE DB KJJ, MJJ, BJJ, USET, GOA, GOAT, KAA, DYNAMICS, KGG , 
PCDBS , AGG, PVT, GPLS , SILS, SLT, DM, CASES, EQEXINS, 
MATPOOL, GM, KOO, LOO, DIT, BGPDTS , CSTMS, EST, MPTS , 
XYCDBS , LLL, KLL, ULL, EMAP, MAPS, ETT, PVTS, KFS , KSS, 
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KSF, KW, EDT, GE0M2S , GE0M3S, ECTS , EPTS , INDTA, KELM, 
EDICT, GPECT, VELEM, EQEXINX, ECTX, BGPDTX , SILX, 

GE0M4S , VGFS , SPECSEL, GPDTS , MELM, MDICT, GE0M1 , 

GEOM2 , GEOM3 , GEOM4 , EPT , MPT, SLIST, VACOMPR, 

$ DYNAMICS ONLY 

BGG, K4GG, MFF , MLR, MLL , MRR, CMLAMA, CMPHA , CMPHO, 
MLAA1 , MR, GOAQ, KLAA , MLAA, MAA, MAR, MEA, BAA, K4AA , 
MGG, 

$ PHASE 2 

PHA, LAMA, PHFA, LAMAF , PHSA, LAMAS, 

$ PHASE 3 

PHG, PHQG, PUG, GPSF , OPHG1 , OPHQG1 , OPHEF1, OPHES1 , 
OPHEE1 , OPHESM1 , OPHESG1 , OPHEEM1 , OPHEEG1 , PHGPS , 
OPHNRGY , OPHGPFB, CASEDR, XYCDBDR, OPHGPS1 , 

$ EQUIV SCRATCH 

MKAA, MMAA , OLB, PCDB, VAFS, XYCDB, EDOMX , CASECC $ 
ALTER 36,36$ 

CALL MODERS MR, USET, DM, CASES, DYNAMICS, MMAA, MKAA, GPLS, 
SILS, EED, EQEXINS , VACOMPR, DMI , DMINDX , DTI, DTINDX/ 
PHA , LAMA/ 

NORSET/ NOLSET/ 'MODES'/ ASING/ FALSE/ FALSE/ NOQSET $ 
COMPILE MODERS SOUIN=MSCSOU NOLIST NOREF 
ALTER 1,1$ 

SUBDMAP MODERS MR, USET, DM, CASES, DYNAMICS, MMAA, MKAA, 
GPLS, SILS, EED, 

EQEXINS, VACOMPR, DMI, DMINDX, DTI, DTINDX/ 

PHIA , LAMA/ 

NORSET/ NOLSET/ READAPP/ ASING/ CYCLIC/ SECND/ NOQSET $ 
ALTER 84$ 

$ Input user specified response levels, amplification 
$ levels, and listing of modes to be correlated. Input via 
$ bulk data DMI and DTI cards. 

DMI IN DMI , DMINDX/PSD , QUAL , ,,,,,,,/$ 

DTI IN DTI , DTINDX/MODELST , ,,,,,,,,/$ 

TYPE PARM , , I , Y , NMODES$ 

TYPE PARM , , I , Y , NTEST$ 

TYPE PARM, ,I,Y,NMAX$ 

TYPE PARM , , RS , Y , GRAV$ 

TYPE PARM , , I , N , COUNT=l $ 

TYPE PARM, ,I,N,PART1$ 

TYPE PARM, ,I,N,PART2$ 

TYPE PARM , , I , N , COUNT1 $ 

TYPE PARM, ,I,N,COUNT2$ 

TYPE PARM, ,RS,N,PI2=1. 57079633$ 

TYPE PARM , ,CS,N, ALPHA $ 

TYPE PARM, ,RS,N,RADI$ 

TYPE PARM, ,CS,N,MRI$ 

FILE CP=OVRWRT /RP=OVRWRT/A 1 1 =0 VRWRT/A 1 2 =OVRWRT $ 

FILE A21 =OVRWRT /A22 =OVRWRT/B 1 1 =0VRWRT/B1 2=OVRWRT $ 

FILE B 2 1 =OVRWRT /B 2 2 =OVRWRT/C 1 1 =0VRWRT/C1 2=OVRWRT $ 

FILE C21=OVRWRT /C2 2=0VRWRT/D1 1=0VRWRT/D1 2=OVRWRT$ 

FILE D2 l=OVRWRT /D2 2=OVRWRT/KINV=OVRWRT/DEF=OVRWRT$ 
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FILE ACCEL=OVRWRT /ACCELI =OVRWRT/MODE I =OVRWRT $ 

FILE DEFT=OVRWRT/XMODE=OVRWRT/COC=OVrWRT$ 

$ Echo user input of response levels, amplification levels, 
$ and modes to be correlated. 

MATPRN PSD , QUAL//$ 

TABPT MODELST//$ 

$ Retrieve and equivalence analysis matrices. GDR 
$ stiffness, GDR mass, eigenvectors, and eigenvalues. 

EQUIVX KXX/KANL/ ALWA Y S $ 

EQUIVX MXX/MANL/ ALWA Y S $ 

EQUIVX PH I X/PH I ANL/ ALWA Y S $ 

LAMX , , LAMA/LAMAT/- 1 $ 

MATPRN KANL , MANL , PHI ANL , LAMAT/ / $ 

$ Partition eigenvector and eigenvalue matrices for the 
$ modes to be correlated. Modes to be correlated are 
$ identified by DTI entry MODELST. There may be at most 
$ as many modes requested as the number specified by 
$ NMODES . 

EQUIVX LAMAT/LAMTMP/ ALWAY S $ 

MATGEN ,/LAMDA/7/l/5$ 

MATGEN ,/PHIRED/7/NTEST/l$ 

C0UNT=1$ 

DO WHILE (COUNT <= NMODES )$ 

PARAML MODELST// ' DTI 7 /COUNT/1 //S , N , PART$ 

PART 1 =PART- 1 $ 

PART 2 =NMAX-PART $ 

IF (PART = 1) THEN$ 

MATGEN ,/CP/6/NMAX/PART/PART2$ 

MATGEN ,/RP/6/NMAX/PART/PART2$ 

ELSE$ 

MATGEN ,/CP/6/NMAX/0/PARTl/l/PART2$ 

MATGEN ,/RP/6/NMAX/0/PARTl/l/PART2$ 

ENDIF$ 

PARTN LAMTMP, ,RP/LAMTP, A21 , A12 , A22/l$ 

PARTN PHIX , CP , /PHITP , B2 1 , B12 , B22/l$ 

IF (COUNT > 1) THEN$ 

COUNTl=COUNT-l$ 

TRNSP LAMDA/TEMP1 $ 

TRNSP LAMTP/TEMP2 $ 

APPEND TEMPI , TEMP 2 /TEMP 3 $ 

TRNSP TEMP 3 /LAMD ATX $ 

EQUIVX LAMD ATX /LAMDA/ALW A Y S $ 

APPEND PHIRED , PHITP/PHITX$ 

EQUIVX PHITX/PHIRED/ALWAYS$ 

ELSE$ 

EQUIVX LAMTP /LAMDA/ ALWAY S $ 

EQUIVX PHITP/PHIRED/ALWAYS$ 

ENDIF$ 

COUNT=COUNT+ 1 $ 

ENDDO$ 

MATPRN LAMDA, PHIRED//$ 

$ Create the nodal directional matrix NDIR by extracting 
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$ each term in PHIRED and selecting its directional 
$ orientation (e.g. +1.0 or -1.0). Directional terms are 
$ assembled in NDIR matrix of size [NTEST X NMODES]. 

EQUIVX PHIRED/NDIR/ ALWAYS? 

COUNT1=0$ 

DO WHILE ( COUNT 1 < NMODES)? 

MATGEN , /CP/6/NMODES/COUNT1/1? 

COUNT2=0? 

DO WHILE (COUNT 2 < NTEST)? 

MATGEN , /RP/6/NTEST/COUNT2/1? 

PARTN NDIR , CP , RP/A1 1 , A2 1 , A1 2 , AOLD/1 ? 

NORM AOLD/ANEW/? 

MERGE All , A21 , A12 , ANEW , CP , RP/NDIR1/1 ? 

EQUIVX NDIR1 /NDIR/ALWAYS? 

COUNT 2 =COUNT 2 + 1 ? 

ENDDO? 

COUNTl=COUNTl+l? 

ENDDO? 

MATPRN NDIR//? 

? Create the nodal acceleration vector ACCEL for each mode 
? by extracting each term from NDIR, PSD, QUAL and LAMDA 
? then solving Miles equation. 

MATGEN , /ACCEL/7/NTEST/NMODES? 

COUNT1=0? 

DO WHILE ( COUNT1 < NMODES)? 

MATGEN , /CP/6/NMODES/COUNT1/1? 

MATGEN , /CPL/ 6/5/2/1/2? 

MATGEN , /RPL/6/NMODES/COUNT1/1? 

PARTN LAMDA, CPL, RPL/, , ,FREQI/1? 

COUNT 2=0? 

DO WHILE ( COUNT 2 < NTEST)? 

MATGEN , /RP/6/NTEST/COUNT2/1? 

PARTN NDIR , CP , RP/A1 1 , A1 2 , A2 1 , NDIRI/1? 

PARTN PSD , CP ,RP/B11 ,B12 ,B21 ,PSDI/1? 

PARTN QUAL , CP , RP/C11 , Cl 2 , C21 , QUALI/1? 

PARTN ACCEL, CP, RP/D11 ,D12 ,D21,DOLD/l? 

SMPYAD QUALI , PSDI , FREQI , , , /RAD/3/1? 

PARAML RAD// # NULL V///S,N,NOPG? 

IF ( NOPG = -1) THEN? 

RADI=0 . 0 
ELSE? 

PARAML RAD//'DMI '/1/1/S, N, RADI? 

ENDIF? 

RADI=( 3 . 0*GRAV*SQRT(PI2*RADI ) ) 

ALPHA=CMPLX ( RADI ) 

ADD NDIRI , /ACCELI/ALPHA? 

EQUIVX ACCELI /DNEW / ALWAY S ? 

MERGE Dll , D12 , D21 , DNEW , CP , RP/ACCEL1/1? 

EQUIVX ACCELI /ACCEL/ ALWAYS? 

COUNT 2 =C0UNT 2 + 1 ? 

ENDDO? 

COUNT 1 =COUNT 1 + 1 ? 
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ENDDO$ 

MATPRN ACCEL//$ 

$ Perform the static solution [KANL] A -1* [MANL] * {ACCEL} for 
$ each mode selected to generate the 'psuedo' static mode 
$ shapes. Each mode shape is mass normalized for the 
$ correlation process. After each mode shape is generated 
$ its influence in the remaining static solution is removed 
$ by performing a deflation of the static matrix in the 
$ form: [ DEF ] = [ KANL ] A - 1 * [ MANL ] - (LAMDA) A -1 *[PHITST] A T 
$ *[PHITST]* [MANL]. Intermediate modes not selected for 
$ correlation are also swept from the system using the 
$ corresponding mode shape provided by analysis. 

MATGEN ,/KINV/7/NTEST/NM0DES? 

MATGEN , /PHITST/7 /NTEST/1$ 

MATGEN ,/M0DEI/7/NTEST/l? 

MATGEN , /DEF/7 /NTEST/NMODES$ 

MATGEN ,/DEFT/7/NTEST/NM0DES? 

SOLVE KANL , /KINV /3 $ 

MPYAD KINV, MANL / /DEF/0/l$ 

C0UNT=1$ 

C0UNT1=1$ 

PARAML MODELST// ' DTI ' /C0UNT1/1//S , N , PART? 

DO WHILE (COUNT <= NMAX )$ 

IF (COUNT = PART) THEN? 

MATMOD ACCEL , , , , , /ACCELI , /1/COUNT1? 

MPYAD DEF , ACCELI ,/XMODE/O/l? 

SMPYAD XMODE , MANL , XMODE , , , /MFSQR/3/1///1/0/0? 

PARAML MFSQR//'DMI ' /1/1/S, N, MR? 

MR=1 . 0/SQRT(MR) ? 

MRI =CMPLX ( MR ) ? 

ADD XMODE, /XTEMP/MRI? 

MESSAGE / / ' CALCULATION OF STATIC MODE SHAPE = '/PART? 

IF (COUNT = 1) THEN? 

EQUIVX XTEMP/PHITST/ALWAYS? 

ELSE? 

APPEND PHITST , XTEMP/XMODE? 

EQUIVX XMODE/PHITST/ALWAYS? 

ENDIF? 

MATMOD PHIANL, , , , ,/MODEI , /1/COUNT? 

SMPYAD MODEI ,MODEI ,MANL, , , /NUMER/3////0/1/0? 

PARAML LAMAT// ' DMI ' /1/COUNT/S , N , EIGVAL? 

MESSAGE //'EIGENVALUE = '/EIGVAL? 

EIGVAL=-1 . 0 /EIGVAL? 

ALPHA=CMPLX ( EIGVAL ) ? 

ADD DEF , NUMER/DEFT/ ( 1 . 0 , 0 . 0 ) /ALPHA/O? 

EQUIVX DEFT/DEF/ALWAYS? 

COUNT=COUNT+l? 

COUNT 1 =COUNT 1 + 1 ? 

PARAML MODELST// ' DTI ' /COUNT1/1//S , N , PART? 

ELSE? 

MATMOD PHIANL , , , , , /MODEI , /1/COUNT? 

SMPYAD MODEI , MODEI , MANL, , , /NUMER/3/1///0/1/0? 
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PARAML LAMAT// ' DMI ' /I /COUNT/S , N , EIGVAL$ 

EIGVAL=-1 . 0/EIGVAL$ 

ALPHA=CMPLX ( EIGVAL ) $ 

ADD DEF , NUMER/DEFT/ ( 1 . 0 , 0 . 0 ) /ALPHA/0$ 

EQUIVX DEFT/DEF/ALWAYS$ 

COUNT=COUNT+ 1 $ 

ENDIF$ 

ENDDO$ 

MATPRN PHITST//$ 

$ Perform the Cross Orthogonality computation for the modes 
$ selected. Orthogonality is provided by [ PHIRED ] A T* 

$ [ MANL ] * [ PHITST ] . 

MATGEN , /COC/7/NTEST/NMODES$ 

SMPYAD PHIRED, MANL, PHITST, , , /COC/3/l///l/0/0$ 

MATPRN COC//$ 

MESSAGE //'CROSS ORTHOGONALITY COMPUTATION COMPLETE '$ 

END$ 



APPENDIX B 


CANTILEVERED BEAM MODEL 

Eigenvalue Analysis, Input 
Eigenvalue Analysis, Output (Abridged) 
Random Response Analysis, Input (Abridged) 
Random Response Analysis, Output (Abridged) 
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SOL 103 
TIME 4 
CEND 

SEALL = ALL 
SOPER = ALL 
SET 100 = 21,41,61 

TITLE = High Fidelity Cantilevered Bean Model 
ECHO = SORT 
SPC = 5 
SOBCASE 1 
METHOD = 1 
BEGIN BOLK 
PARAM AOTOSPC YES 
PARAM COOPMASS 1 

PARAM WTHASS 1. 

PARAM GRDPNT 0 


EIGR 

+ 

PBAR 

CBAR 

1 

A MASS 
1 
1 

SOT 

1 

1 

0 . 

6.-5 

1 

2000. 

5.-10 

2 

1 . 

20 + 
1 . 

0 

CBAR 

2 

1 

2 

3 

1 . 

1 . 

0 

CBAR 

3 

1 

3 

4 
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APPENDIX C 


STDCE EXPERIMENT PACKAGE 

Eigenvalue Analysis, Input (Abridged) 
Eigenvalue Analysis, Output (Abridged) 
Correlation Analysis, Input (Abridged) 
Correlation Analysis, Output (Abridged) 


110 



Ill 


SOL 103 
TIHE 10 
CEND 

TITLE = STDCE EXPERIMENT PACKAGE - AS BCILT FLIGHT HARDWARE 
SUBTITLE = NORMAL MODES ANALYSIS - STATIC CONDESATION EHPLOYED 
ECHO = SORT 
SPC = 41 

VECTOR (PLOT) = ALL 
HETHOD=l 
DYNRED-2 
BEGIN BOLK 

EIGR 1 SINV 0 . 200. 100 + A 

+ A HASS 
DYNRED 2 200. 

$ DOF ASSOCIATED WITH ACCELEROMETER LOCATIONS/DOF 
ASET1 123 603 744 536 149 

ASET1 23 224 

ASET1 2 159 12 

$ PARAMETERS ASSOCIATED WITH EFFECTIVE HASS CALCULATIONS 
PARAM KEPRT 1 

PARAM EFWGT 2 

PARAM CHKSTIF 1 

PARAM CHKMASS 1 

$ 

$ BEGINNING OF EXPERIMENT PACKAGE STRUCTURE (C200) 

$ 

WWWMMMtWMfittttiWtitiMMWWMmttMW 

# 

t DUE TO ITS VOLUMINOUS NATURE, THE BULK DATA INPUT FOR THE 
f EXPERIMENT PACKAGE MODEL HAS BEEN OHITTED. THE COHPLETE MODEL 
f CAN BE OBTAINED FROM THE AUTHOR. 

IwWWMWUMWWmMiWWItmWWMtitWWtt 
$ PARAMETERS ASSOCIATED WITH THE MODAL SOLUTION 
PARAM AUTOSPC YES 

PARAM COUPMASS 1 

PARAM GRDPNT 0 

PARAM K6ROT 1000. 

PARAM WTMASS 1. 

ENDDATA 
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MATRIX EFWFRRCT IS THE FRACTION OF THE EFFECTIVE WEIGHT FOR EACH MODE 
(ONCE AGAIN, FOR THE ASSEMBLY) 



1] 1.0572D-02 -1.5857D-01 1.9509D-01 2.9864D-01 6.1923D-03 -4.6712D-05 

COLUMN 2 ROWS 1 THRU 6 


115 




116 


SOL 103 
TIME 10 

wwwmwtMwwtwtitwwMUMtmw 

CORRELATION DHAP ALTER IS INSERTED HERE 

tWMMWMWfWmmMtWMWMtWWi 

cm 

SET 100 = 603, 744, 536, 224, 149, 159, 12 
TITLE = STDCE EXPERIMENT PACKAGE - AS BUILT FLIGHT HARDWARE 
SUBTITLE = NORMAL HODES ANALYSIS WITH CROSS ORTHOGONALITY CHECK 
SEALL = ALL 
SUPER = ALL 
ECHO = SORT 
SPC = 41 
VECTOR = 100 
METHODS 
DYNRED=2 
BEGIN BULK 

WWMMiWmWltmWWWMMMWWi 

REMAINDER OF INPUT REMAINS UNCHANGED BUT 

HAS BEEN OMITTED HERE FOR BREVITY 

tutttmimmmmtmtmmmmmmitmim 


DMI 

PSD 

0 

2 

1 

0 


16 

5 

DMI 

PSD 

1 

1 

.0301 

3.0142 

.0010 

.0165 

.0002 


.0191 

.8254 

.0010 

.1957 

.0049 

.1695 

.1540 

.0012 

.0083 

.4217 

.0009 

DMI 

PSD 

2 

1 

.0147 

.0013 

.0021 

.0196 

.2488 


.0249 

.0442 

.0866 

.0287 

.0054 

.4217 

.0590 

.4642 

.0049 

.1695 

.5623 

DMI 

PSD 

3 

1 

.0169 

.2738 

.0026 

.0020 

.0178 


.0249 

.0402 

.0909 

.0301 

.0065 

.4870 

.0464 

.5109 

.0054 

.1695 

.6190 

DMI 

PSD 

4 

1 

.0562 

.0649 

.1287 

.0316 

.0001 


.0029 

.1540 

.0075 

.2260 

.9085 

.1334 

1.9573 

.0075 

1.1007 

.3318 

.0506 

DMI 

PSD 

5 

1 

.0162 

.0442 

.0079 

.0000 

.0001 


.0037 

.2054 

.0044 

.0301 

.1957 

.0422 

.0909 

.0054 

.3625 

.0590 

.0037 

DMI 

QUAL 

0 

2 

1 

0 


16 

5 

DMI 

OUAL 

1 

1 

1.9545 

195.727 

.0649 

1.0714 

.0130 


1.2403 .0065 
53.5974 12.7078 

.3182 

11.0065 

10.0000 

.0779 

.5390 

27.3831 

.0584 

DMI 

QUAL 

2 

1 

1.5474 

.1368 

.2211 

2.0632 

26.1895 


2.6211 

4.6526 

9.1158 

3.0211 

.5684 

44.3895 

6.2105 

48.8632 

.5158 

17.8421 59.1895 

MI 

QUAL 

3 

1 

1.7789 

28.8211 

.2737 

.2105 

1.8737 


2.6211 

4.2316 

9.5684 

3.1684 

.6842 

51.2632 

4.8842 

53.7789 

.5684 

17.8421 65.1579 

DMI 

QUAL 

4 

1 

2.9119 

3.3627 

6.6684 

1.6373 

.0001 
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.1503 .3886 47.0725 101.415 .3886 57.0311 17.1917 2.6218 

7.9793 11.7098 6.9119 

DKI QOAL 5 1 .7902 2.1561 .3854 .0000 .0046 

.1805 .2146 9.5463 4.4341 .2634 17.6829 2.8780 .1805 

10.0195 1.4683 2.0585 

DTI,HODELST,0 
DTI , MODELST ,1,1 
DTI,MODELST,2,3 
DT^HODELST^^ 

DTI ,H0DELST / 4 , 5 

DTI,HODELST,5,7 

PARAM,NHODES,5 

PARAM,NTEST,16 

PARAM,NHAX,8 

PARAM,GRAV,386.4 

PARAM HEWSEQ -1 

ENDDATA 
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1 4.9000D-03 1.6950D-01 5.6230D-01 4.4200D-02 2.8700D-02 4.2170D-01 
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11) 5.7031D+01 1.7192D+01 2.6218D+00 7.9793D+00 1.1710D+01 6.9119D+00 

COLUHH 5 ROWS 1 THRU 16 



1) 7.9020D-01 2.1561D+00 3.8540D-01 O.OOOOD+OO 4.6000D-03 1.8050D-01 2.1460D-01 9.54630+00 4.4341D+00 2.6340D-01 
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1) 3.8671D+03 -2.0657D+01 -7 . 39960+01 -7.4112D+00 5.1865D+00 1.9035D+01 -2.3374D+01 3.1034D+01 -5.9529D+01 -4.5280D+02 

11) 3.3613D+02 4.2935D+01 3.4230D+02 8.7708D+00 9.7407D+01 2.18910+01 

COLUMN 2 ROWS 1 THRU 16 — - 



3.5291D+02 1.6954D+04 2.8885D+04 -3.4433D+04 8.4153D+02 -3.5928D+04 





1) 3.3613D+02 2.4570D+03 -1.9O72D+02 6.0953D+02 -7.6812D+01 -5.2188D+02 5.1488D+02 -1.1210D+06 1.6305D+04 4.9368D+05 

11) 1.2568D+06 -3.5982D+04 -4.5721D+05 -1.6238D+04 -3.7406D+04 -3.7624D+04 

COLUMN 12 RCWS 1 THRU 16 
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